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ABSTRACT 

The rapid rotation of Be stars may be caused in some cases by past mass 
and angular momentum accretion in an interacting binary in which the mass 
donor is currently viewed as a small, hot subdwarf stripped of its outer envelope. 
Here we report on the spectroscopic detection of such a subdwarf in the Be 
binary system FY Canis Majoris from the analysis of data acquired by the IUE 
spacecraft and KPNO Coude Feed Telescope over the course of 16 and 21 years, 
respectively. We present a double-lined spectroscopic orbit for the binary based 
upon radial velocities from the IUE spectra and use the orbital solutions with 
a Doppler tomography algorithm to reconstruct the components' UV spectra. 
The subdwarf is hot (T e g = 45 ± 5 kK) and has a mass of about 1.3M & and a 
radius of about 0.6 R Q . It contributes about 4% as much flux as the Be star does 
in the FUV. We also present observations of the Ha and He I A6678 emission 
features that are formed in the circumstellar disk of the Be star. Orbital flux and 
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velocity variations in the He I A6678 profile indicate that much of the emission 
forms along the disk rim facing the hot sub dwarf where the disk is probably 
heated by the incident radiation from the subdwarf. A study of the FUV infall 
shell lines discovered in the 1980s confirms their episodic presence but reveals 
that they tend to be found around both quadrature phases, unlike the pattern 
in Algol binaries. Phase-dependent variations in the UV N V doublet suggest 
the presence of a N-enhanced wind from the subdwarf and a possible shock- 
interaction region between the stars where the subdwarf 's wind collides with the 
disk of the Be star. 

Subject headings: stars: emission-line, Be — stars: individual (HR 2855, HD 58978, 
FY CMa) — stars: binaries: spectroscopic — stars: evolution — stars: subdwarfs 



Introduction 



The component stars of massive interacting binaries experience profound transforma- 
tions over their lifetimes. Their evolutionary stages are define d by how mass and angular 
mom entum are shared and lost during the mass transfer stages (IWellstein. Langer. fc Braun 
200ll ). The originally more massive star grows faster and upon reaching a size comparable 
to its Roche lobe begins mass and angular momentum transfer to the mass gainer compan- 
ion. The orbit shrinks until the mass ratio reverses, but mass transfer will continue with 
an expanding orbit as long as the donor grows in proportion. The result will be a donor 
that is strippe d of most of its out er envelope and a gainer that now appears rejuvenated by 
mass transfer (IDray fc Toutl 120071 ). Mass transfer may also cause the gainer to spin up to a 
rotational speed close to t he critical value where centripetal and gr a vitational a c celera tions 
are equal at the equator ( Petrovic. Langer. &: van der Huchtll2005l ). iPols et al.l (Il99ll ) and 
van Bever &: Vanbeverenl (119971 ) suggest that these mass gainer stars may include some of 
the rapidly rotating, classical Be stars, B-type stars th at display emission lines formed in 
outflowing, circumstellar disks (IPorter &: Riviniud 120031 ) . 



There is now a substantial body of evidence that the supports the idea that some Be 
stars were spun up by mass transfer in interacting binaries. Most of the known massive 
X-ray binar ies consist of a Be star and neutron star companion, the remains of the former 
donor star (jCoell2000l ). Donor stars that end up with a mass lower than the Chandrasekhar 
limit will probably avoid a supernova explosion, but will appear first as a hot subdwarf star 
and then eventually as a massive white dwarf. Although there are no known examples of Be 
plus white dwarf binaries, there are a few identified cases of Be plus hot subdwarf binaries. 
The first discovery of the spectrum of a subdwarf companion was made from the analysis 
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ofUV spectra from the International Ultraviolet Explo rer (IUE) spacecraft (IThaller et al. 
19951 ) and the Hubble Space Telescope ( IGies et al.lll998l ) of the well-known Be star <p Persei. 
The hot companion in <fi Per is an object of solar mass and radius but with an effective 
temperature of 53 kK. However, the hot subdwarf is much fainter than the Be star, and 
its presence is perhaps best observed through i ts heating effects on t h e disk gas nea r est to 
the subdwarf Jstefl. Hummel, fc Rivinkislboool ; iHummel k StefllboOlh . iMaintz et all (bo05h 
show that the binary 59 Cygni shows many of the same characteristics found in Per and is 
probably the second example of a Be plus hot s ubdwarf system. A t hird candidate, HR 2142 
( IPetersI Il983l ). may also belong this category (IWaters et al.lll99ll ) but at the present the 
evidence is not compelling. The fact that the nature of the companions in cf) Per and 59 Cyg 
were only discovered after a century of spectroscopic observation, however, highlights the 
difficulties of their detection (IGiea |2000| ) . 



Here we report on the identification of a hot subdwarf companion in the Be binary 
FY CMa (HR 2855, HP 58978). Balmer emissi on lines in the Be star prim ary in this system 
(BO. 5 I Vnpe. Iffiltner. Garrison, Schildl [l969l: B0.5 IVe. IslettebaS Il982h were first discov- 
ered by lPickering Fleming (119051 ). Although at first glance it appears to be a classical Be 
star, it has historically displayed muc h more spectroscopic variabi l ity in its optical and UV 



spectrum than is typical for Be stars (jSlettebak!ll982l ; iPetersI Il982l . Il988l ; iGrady et al.lll988 



Caoll200ll ). Some of the first IUE observations of FY CMa showed evidence o f sharp, red- 
shifted, " shell" spect ral lines that varied significantly on a timescale of days (IGrady et al. 



19881 ). IPetersI (119881 ) argued that the time scales of variability indicated a binary origin. 



Spectr oscopic variability with a period of 37.26 d was subsequently found by lRivinius et al. 



( 120041 ) who s howed that the streng th and profile of the N V A 1240 wind line fluctuates with 
this period. iRivinius et al.l (120041 ) also presented observations of the Ha and He I A6678 
emission lines that suggested that they also varied in radial velocity and strength in much 
the same way as observed in Per and 59 Cyg, and sugge s ted th at FY CMa might be gener- 
ically similar to these systems. However, IRivinius et al.l (120041 ) could not confirm binarity 
as the orbital elements and the nature of the companion remained unknown. 

Below we present an investigation of FY CMa that is based upon a re-examination of 
the IUE spectra and an analysis of a large collection of optical spectra obtained between 
1985 and 2006 (§2). We first discuss the use of cross-correlation measurements of the IUE 
spectra to obtain radial velocities for both components and to derive the orbital elements 
(§3). We then use the orbital velocity curves to reconstruct the individual UV spectra of 
the components using a Doppler tomography algorithm, and we compare the reconstructed 
spectra to models to estimate basic stellar properties (§4). We then describe the orbital and 
temporal variations we find in the Ha and He I A6678 emission lines (§5). We argue in §6 
that the emission line variations are probably due to heating of the part of the disk rim 
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that faces the hot subdwarf companion. We also discuss the long-term and phase-dependent 
variability in the strength and velocity of the Fe III shell-type infall lines that were discovered 
in the mid-1980s and suggest two possible mechanisms for their formation. Finally, we show 
that phase-dependent profile, strength, and velocity variability in the N V wind lines suggest 
that the subdwarf 's wind is N-enhanced and that there is a significant interaction between 
the subdwarf 's wind and the disk of the Be star. 



2. Observations 

Ninety-seven high resolution FUV spectra were obtained with IUE over the period from 
1979 to 1995. The spectra were secured with the Short Wavelength Prime (SWP) camera 
in echelle mode, and the extracted and flux calibrated spectra from the final archive were 
downloaded from the Multimission Archive at the Space Telescope Science Institute^. We 
binned each spectrum onto a uniform log A wavelength grid (bin size equivalent to 10 km s _1 ), 
aligned them with a global average using cross-correlation shifts from the vicinity of strong 
interstellar lines (after which the interstellar lines are removed by interpolation), and then 
rectified the fluxes by a spline fit to a series of pseudo-continuum zones. We discuss the radial 
velocity measurements and the appearance of the UV spectra in the next two sections. 

We also obtained 158 red spectra of FY CMa with the Kitt Peak National Observa- 
tory 0.9 m Coude Feed Telescope between 1985 and 2006. The various observing runs are 
summarized in Table 1 that lists the heliocentric Julian dates of observation, Besselian year 
range, wavelength range, spectral resolving power, number of individual spectra obtained, 
and a summary of the instrumental configuration. The spectra all record the Ha A6563 
emission line and all but six also record the He I A6678 feature. All the spectra were made 
with the long collimator, camera 5, and either grating A (632 grooves mm -1 ) or grating B 
(312 grooves mm -1 ) in second order. The CCD detector of choice changed over the years 
as better instruments became available, and these are listed in the final column using their 
standard KPNO designations. Most of the spectra are well exposed and have a S/N of 200 
or better in the continuum. 

The spectra were extracted and calibrated using standard routines in IRAF@. All the 
spectra were rectified to a unit continuum by fitting line-free regions with a low-order poly- 



4 http:/ /archive. stsci.edu/iue/ 

5 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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nomial. During most of the runs from 2000 to 2006, we also observed a number of rapidly 
rotating A-type stars (primarily £ Aql) for the removal of atmospheric telluric lines from the 
target spectra. We created a library of atmospheric spectra for each run by removing the 
broad stellar features, and then we divided each target spectrum by the modified atmospheric 
spectrum that most closely matched the target spectrum in a selected region dominated by 
atmospheric absorptions. The spectra from runs before 2000 were transformed to the wave- 
length grid for the fourth to last run listed in Table 1 and then were corrected for telluric 
absorption using the atmospheric spectra from that run (sometimes with less than optimal 
results because of differences in spectral resolution). The spectra from each run were then 
transformed to a common heliocentric wavelength grid. Measurements of these red spectra 
are presented in §5. 



Radial Velocities and Orbital Elements 



The first detection of the UV spectrum of the hot subdwarf in (j) Per was based upon IUE 
spectra t hat recorded many l ines of high excitation transitions especially in the 1310 - 1385 
A range ( Thaller et al.lll995 ). Thus, we began our search for evidence of a hot companion 
by calculating cross- correlation functions (ccfs) of each spectrum with a template spectrum 
within this restricted wavelength range (where the lines of the primary are relatively broad 
and weak). We formed the template spectru m from the synt h etic s pectral atlas associated 
with the non-LTE, line-blanketed models of lLanz fc Hubenyl (120031 ) that were rebinned to 
the IUE wavelength grid. We selected a solar abundance model with T c g = 50 kK and 
log g = 4.5 for the template (close to the optimal values; §4). The resulting ccfs immediately 
revealed a narrow but weak si gnal that displayed a coherent Doppler shift when plotted 
against the orbital period from iRivinius et al.l (120041 ) (see Fig. 1). We then made Gaussian 
fits of these ccfs to determine the radial velocity shifts fo r 50 of the 97 spectra where the ccf 
signal was clearly visible. We used the method of IZuckerl (120031 ) to estimate errors in the ccf 
velocities, and these range from 7 to 20 km s _1 depending on the S/N ratio of the spectrum. 

We also measured ccf radial velocities for the Be star primary using two spectral regions 
where its photospheric lines are particularly strong (1423 - 1523 and 1570 - 1800 A). In 



this case we used a template spectrum derived from the models of lLanz fc Hubenyl (120031 ) 
for T e s = 30 kK and log g = 4.0, first guesses based upon calibrations for the star's spectral 
classification. The resulting ccfs are broad because of the star's large projected rotational 
velocity, and instead of fitting these with a Gaussian, we determined th e radial velocity shifts 
relative to the average ccf. According to the method of IZuckerl (120031 ). the errors in the ccf 
velocities for the Be star are about 3 km s _1 . We then added the Gaussian- fitted radial 
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velocity of the average ccf (V r = 34.3 ±0.3 km s _1 ) to each of these relative velocities to end 
up with the primary star velocities listed in Table 2. This table lists the heliocentric Julian 
and UT dates of mid-exposure, the IUE image number for the spectrum, the orbital phase 
(see below), and the radial velocities V\ and V 2 for the Be star (subscript 1) and the hot 
subdwarf (subscript 2), respectively. Also given with the same subscripts are the observed 
minus calculated residuals (O — C) from the orbital fits (see Table 3), and we see that in 
general the deviations from the fit are comparable to the error estimates and to the scatter 
between measurements of pairs of spectra obtained close in time. 

We made a period search of the two se t s of y elocities using the discrete Fourier trans- 
form met hods of iRoberts. Lehar. fc Dreherl (119871 ). and this analysis confirmed the period 



found by iRivinius et al.l (120041) . We then used the non-linear, least-squares program of 



Morbey &: Brosterhud (119741 ) to determine independent orbital elements for the primary and 
secondary. These produced consistent results for the orbital period and epoch, and we 
adopted the period and epoch from the larger amplitude secondary that had significantly 
lower errors. We found that elliptical solutions did not yield any significant improvement in 
the residuals to the fits, so we assumed that the orbit is circular. Our final orbital elements 
are presented in column 3 of Table 3 and the radial velocity curves are plotted in Figure 2. 
Table 3 lists period, P, epoch of the superior conjunction of the Be star, T S c (defining orbital 
phase 0.0), semiamplitudes, K, and systemic velocities, 7, for both components, mass ratio, 
M2/M1, mass-inclination functions, Msin 3 i, projected semimajor axis, asinz, and r.m.s. 
residuals for each fit (where subscripts 1 and 2 refer to the orbits of the Be and subdwarf 
star, respectively) 



4. Tomographic Reconstruction of UV Spectra 



We used the Doppler tomography algorithm described by iBagnuolo et al.l (119941 ) to 
reconstruct the individual primary and secondary spectra from the composite IUE spectra. 
We took the radial velocity shifts for each component from the orbital solutions in Table 3, 
and then the reconstruction was run for 50 iterations with a gain of 0.8 (the results are 
insensitive to these values). The flux ratio required for the reconstruction was determined 
from the secondary's line depths after making a comparison of the secondary's spectral 
properties with those in model spectra. 

We began our comparis on of the reconstructed secondary spectrum and model spectra 
from lLanz & Hubenyl (120031 ) by considering the appropriate line broadening. We selected 
several strong and minimally blended metallic lines to compare with a grid of synthetic pro- 
files formed by convolution with a rotational broadening function for a linear limb darkening 
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law fjWade & Rucinskil Il985l ) and a range in assumed projected rotational velocity V sini. 
The resulting mean of V sin i = 41 ±5 km s _1 indicates that the secondary's lines are narrow 
and only barely resolved in the IUE spectra. We then experimented with a range of model 
effective temperatures and gravities to try to reproduce the relative patterns of line strength 
(especially in the region surrounding the weak Si III A1299 feature). We also compared the 
reconstructed and model profiles of the Stark broadened He II A1640 feature to help estimate 
the appropriate value of gravity \ogg. The parameters for the best matching solar abundance 
model spectrum are listed in column 3 of Table 4. 

We then estimated the mean UV flux ratio to be F2/F1 = 0.04 ± 0.01 based upon a 
comparison of the line strengths in the best match model and in the reconstructed secondary 
spectrum for a range in test values of flux ratio. We show in Figures 3 and 4 the good 
agreement between the reconstructed secondary and model spectra for two regions of the 
UV spectrum that contain a number of prominent Fe V and O V photospheric lines. 

We also explored the range in synthetic spectral matches that were consistent with 
the appearance of the reconstructed Be star spectrum. Once again we compared the profile 
widths of several conspicuous absorption lines with models for a range in projected rotational 
velocity, and w e found a mean valu e of V sini = 340 ±40 km s~ x . This estimate is similar to 
the result from lFremat et al.l (120051) . 37 5 km s _1 , but it is larger than those reported in other 



optical studies (280 km s" 1 . ISlettebaJd 11982 ; 155 km s' MAbt. Levato. fe Grossol 120021 ) . We 



suspect that the optical lines may at times appear narrower because of unrecognized shell 
components formed in the Be star's disk. We estimated the best match effective temperature 
based on the metallic line strengths (avoiding the use of narrow features like Si III A1299 that 
are probably partially formed in the disk; §6), and the best match model has a temperature 
consistent with the star's spectral classification (Table 4). We simply assumed the gravity 
for the Be star based on its luminosity class since there are no easily made gravity estimates 
from the UV photospheric li nes. We note that both the temperature and the gravity agree 
well with estimates made by iPeterd (119761 ) based upon the continuum colors and the profiles 
of thy and H5. 

FY CMa is not a known eclipsing binary, so we cannot directly determine the system 
inclination and masses. However, the appearance of shel l lines in the UV and optical spectra 
suggests that the inclination is probably large, i > 77° ( ]Hanuschikill996l). if the disk r esides 
in the orbital plane. The typical mass for a B0. 5 V star is ~ 13.2M ( )Harmanedll988l ). and 
from the value of Mi sin 3 i in Table 3, this would imply that the inclination could be as small 
as % = 66°. The range in probable masses is given in Table 4. 

We can estimate the radius of t he Be star primary from its visu al magnitude V and its 

parallax tc as measured by Hipparcos (|ESA||1997I ). Adopting V = 5.63 (jSterken. Vogt. fe Mennickent 
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1996|), E(B - V) = 0.14 (IKaiserl Il989h. Aw = 3.1E(B - V) tt = 2.30 ± 0.70 mas, and a 
bolometric correction of BC = —2.65 (ILanz Hubenyll2003l ). then the radius from the lu- 
minosity and adopted temperature is R± = 5.3i? , where we have ignored the minor flux 
contributions from the disk and subdwarf. The radius of the subdwarf then follows from the 
observed and model UV flux ratios, R2 = 0.63-R©. These mass and radius estimates yield 
gravities that are similar within errors to those from spectroscopic considerations. 



5. Red Spectra Measurements 

The Ha feature in our spectra always appeared as a strong emission line indicating that 
the Be disk was present throughout the period 1985 - 2006. We show in Figure 5 a montage 
of all the Ha profiles arranged according to orbital phase. The profile is often double-peaked, 
usually in the sense that the red peak is stronger around phase = 0.25 while the blue peak 
is stronger nea r phas e = 0.75. Thus, the observations appear to confirm the suggestion by 



Rivinius et al.l (120041 ) that the variations are orbitally modulated. On the other hand, our 
spectra also show that there are temporal variations present unrelated to orbital phase that 
appear as alternating bright and dark horizontal segments in the gray-scale representation. 
We measured the Ha equivalent widths by making a numerical integration of the line flux 
over the interval from 6528 to 6599 A (a region including weak emission components of C II 
A6578,6582), and these are listed in Table 5. The typical equivalent width measurement 
error is 0.04 A. We find that there are significant variations in line strength over time scales 
of a week or so that appear to be unrelated to orbital phase and that probably reflect changes 
in the disk density distribution. 

The outer wings of Ha show a hint of an "S" -shaped velocity variation that is similar 
to the orbital velocity curve of the Be star. Since the outer parts of the profile are formed 
in the disk gas closest to the star, the line wings may share in the orbital motion of the 
underlying star. We measured the bisect or position of the emission line w i ngs ne ar the 25% 



maximum flux level using the method of IShafter. Szkody. fc Thorstensenl ( 119861 ). and these 



wing velocities are listed in column 4 of Table 5 and are plotted together with the UV ccf 
velocities in Figure 2. The typical error for the radial velocities from the Ha wings is 1.0 
km s _1 . The Ha wing velocities indeed follow the orbital motion of the Be star, and the 
orbital elements based upon them are in substantial agreement with those from the UV ccfs 
(see column 2 of Table 2). The differences in systemic velocity 7 are not significant given 
the asymmetries and line blending present in the Ha profile. 



We show the dramatic variations of the He I A6678 profile with orbital phase in Figure 6. 
Here the emission displays a backwards "S" -shaped velocity variation that is similar to the 
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secondary's orbital velocity curve. The emission is strongest around orbital phase 0.0 when 
the hot subdwarf is in the foreground. The emission sometimes assumes a pronounced 
double-peaked shape, especially when the emission is relatively strong. In the same way as 
Ha, the He I line also displays non-orbital related variations that give a banded appearance 
to the gray-scale representation of the variations. We also see the occasional development of 
a central shell absorption feature at phases prior to superior conjunction of the Be star. 

We made simple measurements of the peak velocity positions by locating the maximum 
flux positions in the He I profile. If two peaks were present and clearly separated by a 
local minimum, then we measured both peaks and recorded their mean position V r and 
their separation AV r . If only one peak was visible, then only V r was measured. These 
measurements are given in the final two columns of Table 5 for all those cases where at 
least one strong peak was visible. They are also plotted in Figure 7 together with the 
orbital velocity curve of the hot subdwarf. We see that the velocity amplitude of the He I 
emission is actually larger than that of the subdwarf and that the episodes of greatest peak 
separation occur just prior to the orbital conjunctions. We show in the next section that 
these characteristics are consistent with an emission origin in the outer part of the disk that 
faces the subdwarf. 



6. Discussion 



Our spectroscopic results provide the ingredients to estimate the basic elements of the 
system geometry. Here we discuss the system dimensions and the evidence that the hot 
subdwarf is implicated in local he a ting of th e Be star's disk i n muc h the same way as 
was found for 4> Per by IStefl et al.l (120001 ) and iHummel fc Steffi (120011 ) and for 59 Cyg by 



Maintz et al.l (120051 ) . We present in Figure 8 a sketch of the inferred system geometry as 
viewed from above the orbital plane. For the purpose of this illustration, we assume an 
orbital inclination of i = 70° and a semimajor axis of 112 R Q (where the acceptable range 
is probably 105 — 115i? Q for i = 90° — 66°). We also adopt a Be star radius of Ri = 5.3R & 
(where the probable range is 4.1 — 7.6R Q for distances of 333 — 625 pc). The Be star is 
represented by the solid disk on the left while the hot subdwarf is shown by the small 
dot on the right. The double-lobe surrounding them shows the much larger Roche surface 
boundary. W e can estimate the h a lf ligh t radius of the Ha emitting disk using the method 
described by iGrundstrom fc Gied (120061 ). This estimate is derived from a simple density 
power law for the disk, the mean equivalent width (W\ = —13.6 A), the adopted Be star 
temperature and disk inclination, and an assumed disk outer boundary equal to the volume- 
equivalent, Roche radius. These parameters lead to an estimated disk to star radius ratio 
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of Rd/Ri = 9.4, or a half-light emitting radius of 50R Q for the adopted radius of the Be 
star. We caution that the simplicity of the models may lead to errors in the half-light radius 
of ±30% and that the models assume an azimuthal symmetry for the emission flux, which 
is probably not the case for the disk emission of FY CMa (see below). Nevertheless, this 
simple estim ate of disk radius is ind i cated by the lighter gray shaded region on the left side 



of Figure 8. lOkazaki fc Negueruelal (120011 ) argue that the disks of Be stars in binaries with 
nearly circular orbits are truncated at the 3:1 orbital resonance radius, and this occurs at a 
distance from the Be star of 52R & that is only slightly larger than our disk half-light radius 
estimate. 

We showed in §5 that the He I A6678 emission has an orbital modulation that leads 
to greatest emission strength when the hot subdwarf is in the foreground. Furthermore we 
found that the semiamplitude of the orbital Doppler shift of the emission peaks was larger 
than that of the subdwarf (from Fig. 7, K em « 173 km s _1 > K 2 = 128 km s _1 ). These 
facts suggest that a significant fraction of the emission forms in the part of the disk near 
the subdwarf where the Keplerian velocity of the disk gas is larger than the orbital velocity 
of the subdwarf. Ignoring three-body effects, we can estimate the radius of the emission 
source by assuming both the disk gas and the subdwarf obey Keplerian motion around the 
Be star, R em /a = (M 1 /(M 1 + M 2 )){{K 1 + K 2 )/(K X + K em )f = 0.51 or R em = 58R Q . This 
radius is about 15% larger than the half-light radius, which suggests that the emission forms 
mainly along the outer rim of a disk that may be extended towards the hot subdwarf (see 
Fig. 8). We suspect that this outer disk region is heated by the combined effects of the 
incident subdwarf radiation and the impact of the subdwarf 's wind on the disk. We imagine 
that the rest of the disk also contributes to the line emission flux, and, for example, those 
times when weak emission is observed (see Fig. 7) may correspond to episodes when the gas 
density declines in the outer rim facing the subdwarf and the flux from the rest of the disk 
dominates. 

Since the heated disk gas nearest the hot subdwarf will have a higher Keplerian velocity, 
the center of the heated region will move ahead of the axis joining the stars by an amount that 
depends on the cooling and orbital timescales. We can estimate the phase offset A<f) of this 
asymmetry in two ways. First, the central heated region should reach inferior conjunction 
about midway between the observed velocity extrema at phases <fi = 0.77 and 0.17 (see 
Fig. 7), i.e., at = 0.96. Second, the velocity peak separation is presumably caused by the 
heated region forming in a sector of a ring. We will observe the largest velocity separations 
when the mid-point of the ring segment is oriented along our line of sight. According to the 
peak separations plotted in Figure 7, greatest peak separations occur around = 0.42 and 
cf) = 0.92. This suggests that the apex of the heated rim is offset by about 0.06 in phase (22°) 
ahead of the subdwarf inferior junction at <p = 0.0. From simple geometric arguments, we 
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expect that the maximum half-peak separation V m = max(AK-/2) is related to the emission 
semiamplitude by V m = K em sin 9, where 9 is half the ring sector opening angle. From the 
values of peak separation in Table 5, we estimate V m = 66 km s -1 and 9 = 22°. These 
values were used to plot a representative heated region in the darker gray shaded area of 
Figure 8. The actual extent of the heated ring is probably larger than shown since the 
separated emission peaks actually sample a range of ring azimuthal angles and not just the 
termination angles. 

The subdwarf may also affect the kinematics of the gas motion in the disk rim. We show 
in Figure 9 the orbital phase variations in the appearance of the Si III lines near Si III A 1299. 
The cores of these lines are shell-type features formed in the disk gas seen in projection 
against the photosphere of the Be star. We see evidence that the shell feature broadens from 
= 0.5 to = 0.0 and apparently shifts from a redshifted to a blueshifted position at phases 
surrounding conjunction at = 0.0. We imagine that the motion of the outer disk gas is 
influenced by subdwarf as the disk gas passes and overtakes the star. The gravitational 
attraction would cause acceleration and outward motion as the disk gas approaches the 
subdwarf, but the same forces would cause deceleration and inward motion following closest 
passage. Thus, we would observe infalling (redshifted) disk gas before = 0.0 and then 
outward moving (blueshifted) gas after = 0.0. On the other hand, the disk gas projected 
against the Be star has more nearly Keplerian and tangential motion at = 0.5 leading 
to the appearance of narrow shell lines at that phase. This kind of perturbation in the 
Keplerian motions may also partially explain the appearance of a redshifted shell feature in 
the He I A6678 profiles near = 0.9 (Fig. 6), an orientation where the disk gas projected 
against the Be star would be infalling and relatively dense. 

In addition to the shell features discussed above that predominantly are formed in the 
disk, another circumstellar component is seen in some intermediately-ionized species (e.g. 
S III, Fe III) . These variable features, which rarely show velocities less than the photospheric 
value and can display velocities as large as 50 km s _1 relati ye to the photosphere of the Be 



star, were discovered in IUE images obtained in the 1980s (jPetera Il988l ; iGrady et al.lll988l ) 
and they clearly show the presence of significant infall of material toward the Be star. Often 
times these l ines a re asymmetric with enhanced absorption on the red side of the profile. 



Grady et al.l (119881 ) identified six infall episodes that are best seen in the Fe III (multiplet 34) 
line at 1895.456 A, since the feature does not have to compete with a strong photospheric 
component and it is not blended with interstellar or other circumstellar lines. 

From the entire set of 97 IUE SWP high resolution (HIRES) images that span more 
than 15 years of observation and our knowledge of the orbit of the system, we now have more 
insight on the nature of the infall lines. Their behavior in strength and velocity are shown 
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in Figures 10, 11, and 12. Representative examples of the three basic types of observed 
features, photospheric, photospheric with additional absorption from the Be star's disk, and 
infall are given in Figure 10. In order to eliminate phase-dependent effects, we have chosen 
spectra taken at approximately the same quadrature phase on the trailing hemisphere of 
the Be star but at different epochs. Note that the central core of the profile taken during 
a period of enhanced disk absorption shows the same velocity as the photospheric line, but 
that the feature formed in infalling material has a significant redshift. In the upper two-part 
panel of Figure 11 we plot the residual intensit}{§ of the deepest part of the Fe III line and 
its radial velocity versus observation date. All data are included in the intensity plot, but in 
the velocity plot we use a different symbol for observations that show no, or minimal shell 
absorption. Points that fall above the dashed horizontal line in the upper panel visually 
appeared to be mostly photospheric. The double horizontal dashed lines in the second panel 
delineate the domain of photospheric motion. Features that fall above the double line are 
clearly associated with infall phases. Seven or eight clear infall episodes are observed. The 
data suggest that the duration of the infall phases varies from days to a year, but significant 
data gaps from 1979-85 and 1990-94 render it impossible to determine a firm value or its 
degree of variability. An apparently extended infall phase of about 300 days was seen from 
HJD 2446919-2447215. In the lower two panels of Figure 11, the same data are plotted 
versus orbital phase. From the third panel it can be seen that on the average the shell- 
type absorption is enhanced in the phase interval 0.3-0.7. But from the display in the fourth 
panel it is apparent that features seen around phase 0.5 tend to display the systemic velocity, 
implying they are formed in the disk of the Be star. In the fourth panel core velocities from 
only the spectra showing enhanced shell absorption are compared with the radial velocity 
curve of the Be star. It is clear that infall features with velocities up to 50 km s -1 relative 
to the photosphere are seen at most phases with a minimal occurrence around phase 0.5. 
The infall features tend to cluster around the quadrature points at which their velocities 
are highest relative to the Be star's photosphere. Many of the phase dependent strength 
and velocity variations discussed above can also be seen in the gray-scale plot presented in 
Figure 12. 

The episodic i nfall of mate rial is fundamentally different from that observed in the 



Algol- type binaries (jPeters 1120011 ). where mass motion toward a mass gainer is typically seen 
only between phases 0.75-0.95. In the case of FY CMa, the infall features predominantly 
cluster around the quadrature phases from 0.15-0.35 and 0.65-0.90. The velocities of the 
shell features around phase 0.5 tend to be either photospheric or negative relative to the 



6 The observed flux divided by the flux in the local continuum outside of the spectral feature. 
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Be stafl. iPetera ( 119881 ) discussed various scenarios to explain the observed simultaneous 
presence of infall and outflow features in this system includin g a sudden accretion event, the 
possible presence of magnetic loops, or flow from a polar jet. ICaol (120011 ) reported that infall 
components in the N V doublet were contemporaneously observed with those in He I A5876 
and also suggested that magnetic loops might be involved, but the phase-dependence of the 
infall features does not support this speculation. The infall is more likely a consequence 
of an interaction between the wind from the subdwarf and the periphery of the massive 
disk of the Be star or disk perturbations due to the tidal force produced by the subdwarf. 
The motion of a shock-interaction region would result from variations in the location of the 
interface where there is a balance between the dynamic pressure of the subdwarf 's wind and 
the gas pressure at the periphery of the Be star's disk. It is not clear at this time whether it 
is the subdwarf or the Be star that might be driving the variability though. Epochs when no 
enhanced shell absorption is observed would correspond to minimal wind-disk interaction. 
Alternatively the episodic infall features could simply be a result of gas dynamics in the disk 
due to a tidal effect from the subdwarf. Hydrodynamical sim ulations^] of the gas motion in 
a perturbed disk in a Be/X-Ray binary (jOkazaki et al. Il2002l ) reveal turbulence in the outer 
disk and even the formation of spiral structure. Similar disk variability could be operating in 
FY CMa. The details are probably complex but an interplay of such gravitationally related 
forces could create time variable and sometimes coherent gas flows. For example, the tidal 
deformation of the disk leads to an axial elongation that moves ahead of the motion binary. 
Disk gas returning from the extended regions would be observed to be moving inwards when 
viewed at the quadrature phases. 

Finally, we retu rn to the cyclic variation in the N V A1240 wind line that first led 
Rivinius et al.l (120041 ) to find the orbital period of the binary. The plot of these variations 
in Figure 13 shows that the greatest extent of blueshifted absorption occurs around = 0.0 
when the subdwarf is in the f oreground. We note that such a strong N V wind feature 
is quite unusual for Be stars (IGrady. Bjorkman. fc Snowl Il987l ). and the orbital phasing 
suggests instead that the feature may form through absorption of the Be star's flux by 
the wind of the hot subdwarf when it is in the foreground. The profile would mainly be 
affected by blueshifted wind absorption since any redshifted absorption component would 
probably by canceled by redshifted emission in the subdwarf's wind. This kind of orbital 
wind variation is not seen in the other major wind line, the C IV A1550 doublet, and this 
may suggest that the subdwarf is N-enriched and C-depleted as expected for a star that 



7 Evidence for material outflow is often seen in interacting binaries of the Algol type (jPeters 1120011 ) and 
similar hydrodynamics might be occuring in this Be + subdwarf system. 



3 cf. http: / /harmas. arc. hokkai-s-u.ac.jp/%7Eokazaki/BeX/sim/index. html 
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has been stripped down to near the CNO-burning core. One problem with this explanation 
is that the largest wind speeds observed, ~ 400 km s _1 , are much smaller than the wind 
terminal v elocity of a few thousand km s _1 expected for a radiatively driven wind from the 



subdwarf (ILamers fc Cassinellil Il999l ) . However, it may be that the wind optical depths are 



too small for absorption to occur at higher outflow velocities. 

In 1995 March a sequence of daily IUE SWP HIRES spectra of FY CMa were obtained 
over the course of 16 days. Fortuitously these spectra were centered in time around orbital 
phase 0.0 and extended from phases 0.784 to 0.213. The temporal variations in the N V A1240 
doublet are illustrated in a gray-scale diagram in Figure 14. More detail can be seen from 
the plot of selected observations shown in Figure 15. The N V features initially showed very 
little structure until a core with the velocity of the center-of-mass of the system appeared 
around phase 0.97. The core, still with no velocity shift, reached its greatest strength at 
conjunction. It broadened by phase 0.08, and had completely vanished by phase 0.21. The 
core is probably formed in a shock-heated region between the stars where the N-enhanced 
wind from the subdwarf collides with the disk of the Be star. 

Our work demonstrates that FY CMa represents the third known case of a Be star 
with a hot subdwarf companion, the remains of an intense binary interaction that stripped 
down the mass donor to a tiny fraction of its original mass and left the mass gainer spinning 
rapidly. The subdwarf secondary in FY CMa has a mass that is close to the Chandrasekhar 
limit, and lacking a significant H-envelope, it may represent the progenitor of a supernova 
type SN Ibc. The system is now wide enough that it is unlikely that the Be star will be 
substantially spun down by tidal effects, so its rapid rotation will be important throughout its 
subsequent evolution. Such rapidly rot ating massive stars m ay be related to the progenitors 



of collapsars and gamma-ray bursters (ICantiello et al.l 120071 ). 
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Table 1. Journal of Spectroscopy 



Dates 


Bessclian 


Range 


Resolving Power 


Number of 


Telescope / 


(HJD-2,400,000) 


Year 


(A) 


(A/AA) 


Spectra 


Grating/Detector 


44170 


- 49804 


1979.8 


- 1995.2 


1200 - 


1900 


10000 


97 


/rayEchciic/swp 


46123 


1985.2 


6477 - 


6654 


13000 


1 


KPNO CF/B/TI3 


46153 


1985.2 


6514 - 


6691 


13000 


1 


KPNO CF/B/TI3 


46735 


1986.8 


6511 - 


6688 


13000 


1 


KPNO CF/B/TI3 


46903 


- 46906 


1987.3 


- 1987.3 


6522 - 


6698 


13000 


3 


KPNO CF/B/TI3 


46918 


- 46922 


1987.3 


- 1987.3 


6520 - 


6698 


13000 


6 


KPNO CF/B/TI3 


47033 


- 47034 


1987.6 


- 1987.7 


6518 - 


6695 


13000 


2 


KPNO CF/B/TI3 


47303 


- 47307 


1988.4 


- 1988.4 


6517 - 


6694 


13000 


9 


KPNO CF/B/TI3 


47470 


- 47473 


1988.8 


- 1988.9 


6495 - 


6717 


11000 


11 


KPNO CF/B/RCA2 


47561 


- 47562 


1989.1 


- 1989.1 


6514 - 


6690 


14000 


3 


KPNO CF/B/TI3 


47637 


- 47641 


1989.3 


- 1989.3 


6520 - 


6697 


14000 


4 


KPNO CF/B/TI3 


47940 


- 47941 


1990.1 


- 1990.1 


6523 - 


6699 


14000 


2 


KPNO CF/B/TI3 


47982 


- 47984 


1990.2 


- 1990.3 


6520 - 


6697 


14000 


2 


KPNO CF/B/TI3 


48314 


- 48315 


1991.2 


- 1991.2 


6542 - 


6690 


18000 


3 


KPNO CF/A/ST1K 


48319 


- 48322 


1991.2 


- 1991.2 


6520 - 


6696 


14000 


4 


KPNO CF/B/TI3 


48515 


- 48517 


1991.7 


- 1991.7 


6517 - 


6693 


14000 


5 


KPNO CF/B/TI3 


49058 


- 49061 


1993.2 


- 1993.2 


6384 - 


6723 


19000 


3 


KPNO CF/A/T2KB 


49444 


- 49448 


1994.2 


- 1994.3 


6525 - 


6693 


18000 


4 


KPNO CF/A/T1KA 


51123 


- 51126 


1998.8 


- 1998.9 


6359 - 


6677 


20000 


4 


KPNO CF/A/F3KB 


51192 


- 51197 


1999.0 


- 1999.0 


6420 - 


6737 


20000 


5 


KPNO CF/A/F3KB 


51502 


- 51512 


1999.9 


- 1999.9 


6415 - 


6733 


20000 


16 


KPNO CF/A/F3KB 


51613 


- 51614 


2000.2 


- 2000.2 


6459 - 


6776 


20000 


8 


KPNO CF/A/F3KB 


51615 


- 51618 


2000.2 


- 2000.2 


6413 - 


6730 


20000 


9 


KPNO CF/A/F3KB 


51818 


- 51831 


2000.7 


- 2000.8 


6440 - 


7105 


9500 


9 


KPNO CF/B/F3KB 


51850 


- 51852 


2000.8 


- 2000.8 


6413 - 


6733 


20000 


7 


KPNO CF/A/F3KB 


51890 


- 51902 


2000.9 


- 2001.0 


6440 - 


7105 


9500 


12 


KPNO CF/B/F3KB 


51914 


- 51918 


2001.0 


- 2001.0 


6432 - 


6749 


20000 


18 


KPNO CF/A/F3KB 


53292 


- 53295 


2004.8 


- 2004.8 


6470 - 


7140 


9500 


4 


KPNO CF/B/T2KB 


54021 


- 54025 


2006.8 


- 2006.8 


6473 - 


7143 


9500 


2 


KPNO CF/B/T2KB 



Table 2. IUE Radial Velocity Measurements 



Date UT SWP Orbital Vi (O - C)i V 2 (O - C) 2 

(HJD-2,400,000) (yyyy-mm-dd) Number Phase (km s" 1 ) (km s" 1 ) (km s" 1 ) (km s" 1 ) 



44170.044 1979-10-23 6963 0.985 51.5 16.6 1.6 -17.9 

44330.226 1980-03-31 8617 0.285 23.9 4.4 159.8 3.4 

44853.156 1981-09-05 14910 0.321 22.2 1.6 138.6 -8.4 

44867.935 1981-09-20 15053 0.717 40.5 -7.3 -123.5 -29.2 

44920.747 1981-11-12 15478 0.135 39.0 16.2 



Note. — A machine readable version of the full Table 2 is available in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 
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Table 3. Orbital Elements for FY CMa 



Element 


Ha Wings 


IUE ccfs 


P (days) 


37.253 ± 0.007 


37.257 ±0.003 


T sc (HJD-2,448,000) 


530.19 ±0.76 


529.64 ±0.15 


K x (km s" 1 ) 


17.4 ±0.9 


14.4 ±0.9 


K 2 (km s- 1 ) 




128.2 ±2.2 


71 (km s _1 ) 


47.1 ±0.6 


33.6 ±0.7 


72 (km s _1 ) 




31.2 ±1.7 


M 2 /M 1 




0.113 ±0.007 


Mxsin 3 ^ (M ) 




10.1 ±0.4 


M 2 sin 3 i (M ) 




1.13 ±0.08 


asini (Rq) 




105.0 ± 1.7 


rmsi (km s^ 1 ) 


7.6 


6.5 


rms2 (km s _1 ) 




14.5 



Table 4. Stellar Parameters for FY CMa 



Parameter 


Primary 


Secondary 


T eff (kK) 


27.5 ±3.0 


45 ±5 


log g (cgs) 


4.0 


4.3 ±0.6 


V sini (km s _1 ) 


340 ± 40 


41 ±5 


M/Mq 


10 - 13 


1.1 - 1.5 


Rj Rq 


4.1 - 7.6 


0.5-0.9 
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Table 5. Red Spectra Measurements 



Date Orbital W x (Ha) V(Ka wings) V(Ue I peaks) AV(He I peaks) 

(HJD-2,400,000) Phase (A) (km s" 1 ) (km s" 1 ) (km s" 1 ) 



46122.805 0.399 -10.98 26.8 

46152.702 0.201 -13.32 21.8 216 

46734.931 0.829 -11.22 71.1 -154 

46902.728 0.333 -11.76 30.2 

46903.703 0.359 -12.39 25.2 



Note. — Table 5 is available in its entirely in the electronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 
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Fig. 1. — The orbital phase variations in the cross-correlation functions of the spectra in 
the 1310 — 1385 A region with a hot star template spectrum. The ccfs are shown as linear 
plots (top panel) and as a gray-scale image (lower panel). The intensity in the gray-scale 
image is assigned one of 16 gray levels based on its value between the minimum (dark) and 
maximum (bright) observed values. The intensity between observed ccfs is calculated by a 
linear interpolation between the closest observed phases (shown by arrows along the right 
axis). The weak signal from the hot subdwarf spectrum appears as a backwards "S" feature 
in the lower panel. The scale bar at top right indicates the amplitude of the ccf variation 
relative to the mean sum of the squared differences far away from the optimal correlation 
shift. 
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Fig. 2. — Radial velocity curves for the Be star and its companion. Orbital phase 0.0 corre- 
sponds to the Be star superior conjunction time. The solid circles represent the IUE cross- 
correlation velocities for the Be star while the open circles show the IUE cross-correlation 
results for the hot subdwarf. The plus signs indicate the velocities from measurements of 
the Ha wing bisector positions. The typical radial velocity errors are 3 and 7 km s _1 for 
measurements of the Be star and hot subdwarf, respectively. 
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Fig. 5. — The orbital phase variations in the Ha emission line in the spectra of FY CMa are 
shown in linear plots (top panel) and as a gray-scale image (lower panel). The intensity in the 
gray-scale image is assigned one of 16 gray levels based on its value between the minimum 
(dark) and maximum (bright) observed values. The intensity between observed spectra is 
calculated by a linear interpolation between the closest observed phases (shown by arrows 
along the right axis). The scale bar at top right indicates the spectral flux relative to the 
local continuum flux. 
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Fig. 6. — The orbital phase variations in the He I A6678 emission line plotted in the same 
format as Fig. 5. 



-28- 



300 F 




-200 b , i , , , , i , , , , i 

0.0 0.5 1.0 

ORBITAL PHASE 



Fig. 7. — The radial velocities of the prominent emission peaks in the profile of He I A6678 
plotted as a function of orbital phase. Single dots represent measurements of one peak 
(possibly blended) while dots with vertical segments represent the midpoint velocity where 
two prominent peaks were observed. The end points of these line segments correspond to 
the velocity of each resolved peak. The dotted line represents the orbital velocity curve of 
the hot sub dwarf. 
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Fig. 8. — An empirical model of the system geometry as viewed from above the orbital plane. 
Each tick mark corresponds to a length of 0.5a. The Be primary (left) is surrounded by an 
equatorial disk (lightly shaded region) that is heated along the edge (darker shaded region) 
facing the subdwarf secondary (right). The outer boundary of the Be star disk is probably 
within the critical Roche surface (solid curve). The center-of-mass of the system is indicated 
by the vertical tick mark on the line-of-centers. The orbital phases are indicated on the 
periphery. 
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Fig. 9. — The orbital phase variations in Si III A 1299 and other nearby Si III lines plotted in 
the same format as Fig. 5. These features are a composite of broad photospheric lines and 
narrower "shell" features formed in the disk. 
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Fig. 10. — Selected profiles of the Fe III A1895 feature, all observed near phase 0.75 but at 
different epochs. Examples include an apparently pure, rotationally-broadened photospheric 
profile (dotted line), a profile with weak shell absorption (thin, solid line), and an observation 
made during a prominent infall phase (thick, solid line). The vertical line above the profiles 
shows the observed line center that includes the velocity shift of the Be star. 
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Fig. 11. — The absorption strength and velocity of the core of the Fe III A 1895 A line versus 
time and orbital phase. Panel 1 : The horizontal dashed line represents the core intensity 
above which the feature appears to be primarily photospheric. Panel 2: The double dashed 
lines represent the range in photospheric velocity of the Be star. The observations shown 
with the open diamonds showed either no or minimal shell components. All observations 
above the top line were obtained during infall epochs. Panel 3: The dashed line has the 
same meaning as in panel 1. Enhanced shell absorption is seen between phases 0.3 and 0.7. 
Panel 4 '■ The velocity of the center-of-mass of the system is indicted by a dashed line and 
the radial velocity curve of the Be star with a solid line. Only observations showing a shell 
component are included. Note the obvious presence of infalling material at most phases 
except perhaps near phase 0.5 where material outflow is commonly seen in close binaries 
with early type primaries (Peters 2001). The errors for all observations are less than or 
equal to the dimension of the plotted point. 
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Fig. 12. — The orbital phase variations in Fe III A1895 line plotted in the same format as 
Fig. 5. A narrow "shell" feature can be seen around orbital phase 0.5. 



-34- 




Fig. 13. — The orbital phase variations in the N V A1240 doublet plotted in the same format 
as Fig. 5 (in the velocity frame of the blue member of the doublet). The strong absorption 
occurs when the Be star is viewed through the wind of the hot subdwarf. 
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Fig. 14. — The temporal variations in the N V A 1240 doublet observed in a daily sequence 
of IUE observations obtained in 1995 March. The lower portion shows the average profile 
while the upper panel shows a gray-scale representation of the deviations from the average 
as a function of time. A shell feature appeared near orbital phase 0.0 at HJD 2449796.4. 
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Fig. 15. — Selected profiles of the N V doublet observed during the 1995 March sequence 
show the details of the development and fading of the N V shell component. The feature was 
strongest at inferior conjunction of the subdwarf and essentially absent near the quadrature 
points. The vertical dotted lines show the locations of the observed line centers in the velocity 
frame of the center-of-mass of the system. The y-axis flux offsets are (from top to bottom) 
0, —4, —10, —15, and —20 in the flux units plotted. 
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DVR(HE) 
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Note (1) : HRV(HA) values omitted in Fig. 1 and in the orbital solution of 
Table 3, column 2, for HJD-2400000 = 47639.671 and 47983.611. 
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